Fracturing graphene by chlorination: a theoretical viewpoint 
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Motivated by the recent photochlorination experiment [B. Li et a/., ACS Nano 5, 5957 (2011)], we 
study theoretically the interaction of chlorine with graphene. In previous theoretical studies, cova- 
lent binding between chlorine and carbon atoms has been elusive upon adsorption to the graphene 
basal plane. Interestingly, in their recent experiment, Li et at. interpreted their data in terms of 
chemical bonding of chlorine on top of the graphene plane, associated with a change from sp 2 to 
sp 3 in carbon hybridization and formation of graphene nanodomains. We study the hypothesis that 
these domains are actually fractured graphene with chlorinated edges, and compare the energetics 
of chlorine-containing graphene edge terminations, both in zigzag and armchair directions, to chlo- 
rine adsorption onto infinite graphene. Our results indicate that edge chlorination is favored over 
adsorption in the experimental conditions with radical atomic chlorine and that edge chlorination 
with sp 3 -hybridized edge carbons is stable also in ambient conditions. An ah initio thermodynamical 
analysis shows that the presence of chlorine is able to break the pristine graphene layer. Finally, we 
discuss the possible effects of the silicon dioxide substrate on the chlorination of graphene. 



I. INTRODUCTION 

Graphene, the two-dimensional allotrope of carbon, 
has been a topic of intense study ever since its experi- 
mental discovery in 2004. As the properties of pristine 
graphene are becoming increasingly well understood, the 
research community has turned its attention toward the 
modification of graphene properties. Graphene is a zero 
band gap material that has a linear dispersion in the chi- 
ral electronic bands near the Fermi level, which causes 
some of its extraordinary properties such as high electron 
mobility. The lack of a band gap in bulk graphene poses, 
however, limitations for its use in some applications, such 
as in graphene-based transistors P Cutting graphene into 
narrow ribbons introduces a band gapP An alternative 
way of introducing a gap is chemical modification that 
occurs either at graphene edges, or through adsorption 
or chemical binding to the carbon plane. For instance, 
through hydrogerP^ or fluorine^H^ attachment onto the 
basal plane of graphene a band gap can be induced. Fur- 
thermore, the properties of graphene nanoribbons can be 
modified through edge-terminating group^^D and even 
half-metallicity can enter through inequivalent edgesP^^ 

This far, only the lightest one of the halogen atoms, 
fluorine, has been extensively studied in the context of 
interaction with graphene, both theoreticall}B^3] an( j 
experiment ally!"™ The structures of fully fluorinated 
graphene and hydrogenated graphene are similar; in both 
the covalently bound heteroatoms changing the carbon 
hybridization from sr? to sp 3 and strongly corrugating 
the graphene layer P™ Instead, calculations with heav- 
ier halogens (CI, Br, and I) on graphene with all car- 
bon atoms halogenated have predicted ionic binding to 
graphene with no change in the carbon hybridization, 
retaining the planarity of the carbon network, and no 
induced band gap! 19 | 2Q | The carbon-halogen distances for 
these structures are around 4 A, much higher than typ- 
ical covalent bond lengths. Additionally, studies on ad- 



sorption of both isolated heavier halogen atoms^^H anc [ 
halogen molecules^ on graphene have predicted no co- 
valent bonding between the halogen atoms and graphene 
carbons. For graphene nanoribbons (GNRs), some struc- 
tures with halogen-carbon binding at the edges of a 
zigzag GNR (ZGNR) at low edge coverage have been 
reported.^ 

Recently, Li et alW* studied graphene-chlorine interac- 
tions by using radiation to break CI2 molecules to reac- 
tive chlorine radicals. Based on XPS and Raman spec- 
tra, they report graphene basal plane chlorination as- 
sociated with sp 3 hybridized carbon atoms, the forma- 
tion of graphene nanodomains that are approximately 
30-50 nm in lateral dimension, the appearance of a band 
gap greater than 45 meV, and estimate the degree of chlo- 
rination to 8 atom-%. In another recent experimental 
report, Wu et a/P^ report two different reaction regimes 
for graphene exposed to CI2 plasma. For short exposure 
times, the conductivity of graphene is increased related 
to p-doping, only a weak defect-related Raman peak is 
seen, and pristine graphene is recovered upon annealing. 
For long exposure, an irreversible change in the graphene 
properties is reported, associated with formation of chlo- 
rine binding defect patches as seen in transmission elec- 
tron microscopy images. Thus, in both of these experi- 
ments, changes in the structure of graphene are observed 
upon chlorination. 

The differences in these experiments may be due to 
different chlorination procedure and graphene quality. 
The purpose of this study is to address this discrep- 
ancy between the binding of CI on graphene in the 
photochlorinat ion ex periment 2 ^ and previous theoretical 
calculations l 19 n 21 l 24 l commenting also the very recent ex- 
perimental observations by Wu et al. in Ref. |26| We pro- 
pose that the presence of atomic chlorine can lead to the 
formation of edges in graphene, so that the nanodomains 
observed by Li et a/P^ could actually be chlorine-edged 
nanostructures, and compare the chlorinated edge forma- 



2 



TABLE I. Chlorine binding energies [Eq. |T]), in eV per CI atom] and chlorine equilibrium distance from the carbon plane 
for different supercell sizes and initial positions of the chlorine atom. The calculations starting from the hollow position that 
converge to top positions are excluded from the table. 

bridge hollow top 

Eb,c\ Eb,c\ 2 d Eb,c\ Eb,c\ 2 d Eb,c\ Eb,c\ 2 d 
(eV/Cl) (eV/Cl) (A) (eV/Cl) (eV/Cl) (A) (eV/Cl) (eV/Cl) (A) 



4x4 


-1.040 


0.703 


2.93 








-1.049 


0.695 


2.86 


6x6 


-1.159 


0.584 


3.10 


-1.142 


0.601 


3.23 


-1.161 


0.582 


3.05 


8x8 


-1.210 


0.533 


3.12 


-1.196 


0.547 


3.22 


-1.212 


0.532 


3.08 


12x12 


-1.244 


0.499 


3.17 








-1.248 


0.496 


3.12 



tion to adsorption onto bulk graphene. By considering a 
large number of chlorine-containing edge terminations in 
our density- functional theory calculations, we show that 
chlorine atoms preferably bind to the graphene edges, 
rather than onto the basal plane of graphene. Using an 
ab initio thermodynamics approach, we find that only 
chlorinated armchair (AC) edges are stable with respect 
to pristine graphene, and that in the presence of radical 
chlorine, graphene may spontaneously break into struc- 
tures with chlorinated armchair edges. Additionally, we 
discuss the possible effects that a Si02 substrate can have 
on chlorination. 



II. COMPUTATIONAL DETAILS 

Our calculations were performed using the density- 
functional theory (DFT) with a van der Waals (vdW) 
correction} 2 -^ implemented in the all-electron density- 
functional code, developed at the Fritz-Haber Institute 
(FHI), called FHI-aims. 28 Double numeric plus polariza- 
tion basis set of numerical atom-centered orbital basis 
functions and the PBE 29 exchange correlation function- 
al were used. The structural relaxation (with the vdW 
correction) was converged until the forces acting on the 
atoms were less than 0.001 eV/A and the electronic de- 
grees of freedom were converged below 10 -6 eV. Unless 
otherwise mentioned, no spin polarization was taken into 
account, as in previous calculations chlorine adsorption 
has not been shown to induce spin polarization.^ 

The optimized lattice constant, 1.42 A, of an infinite 
graphene sheet was used in all calculations in the absence 
of a substrate. In order to avoid interaction between the 
periodic images, a vacuum layer of 15 A was placed be- 
tween the graphene layers. The Brillouin zone sampling 
was adjusted to the supercell size, corresponding to a 
48x48 k-point mesh for the graphene primitive unit cell. 
In the calculation of binding energies, the reference state 
for carbon was chosen as the infinite pristine graphene 
sheet, and for hydrogen a hydrogen molecule in the gas 
phase. For chlorine, both molecular and atomic refer- 
ence states were used, as in the photoclorination experi- 
ment radiation is used to break the clorine molecules into 
atomic chlorine radicals. 25 Thus, the formation and bind- 



ing energies calculated with respect to atomic chlorine 
corresponds to reaction conditions, and the molecular to 
ambient conditions. 



III. RESULTS 
A. Chlorine adsorption on infinite graphene 

First, we study the binding of one chlorine atom on the 
infinite graphene sheet. In previous calculations, no cova- 
lent binding between graphene carbon and chlorine atoms 
that could change the hybridization of carbon atoms from 
sp 3 to sp 2 has been found. Instead, the adsorption was 
of ionic type, associated with charge transfer between 
graphene and the adsorbate, that has the lowest-energy 
adsorption positions on top of a carbon atom p^ or on the 
bridge site between two carbon atoms The carbon- 
chlorine distances were around 3 A, 24 clearly longer than 
the covalent bond length between carbon and clorine, 
1.72-1.85 APThe adsorption energies in these calcula- 
tions were -0.3 eV/Cl 24 or -0.80 eV/ClJ^the chlorine 
reference state being atomic chlorine. 

In our calculations, we consider three initial positions 
for chlorine atoms: on top of a carbon atom (t), bridge 
between two carbon atoms (b) and in the center of a car- 
bon hexagon (h). As periodic supercells are used in the 
calculation, we address the finite size effects by calculat- 
ing the binding energy for the single CI atom absorption 
to graphene for different supercell sizes. The binding en- 
ergy is given by 

Eb,C\ = ^Cl-gra - ^gra ~ ^Cb (1) 

where £?ci-gra is the total energy of the graphene su- 
percell with a single chlorine atom, E gra the energy of 
pristine graphene in a supercell of the same size, and Eqi 
the energy of a chlorine atom. When molecular reference 
for chlorine is used, Eqi is substituted by Eq\ 2 /2. 

Table U shows our results for the atomic chlorine bind- 
ing energy and the equilibrium distance for different su- 
percell sizes, ranging from a 4x4 supercell (32 carbon 
atoms) also used in Ref. [3TJ to a 12x12 supercell (288 
carbon atoms). Our results qualitatively agree with those 
of Ref. [211 giving little difference in adsorption energies 
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TABLE II. The chlorination of graphene on S1O2 for dif- 
ferent substrate terminations in the 2x2 and 4x4 graphene 
supercells. The C-Cl bond distance for adsorbed chlorine 
(dc-ci) and the chlorine binding energy with respect to rad- 
ical (Eb, ci) and molecular chlorine (Eb,c\ 2 ) are given for all 
surface terminations. 



termination 


Eb, ci 


Eb, ci 2 


dc-ci 




(eV/Cl) (eV/Cl) 


(A) 


02x2 


-0.920 


0.530 


1.83 


04x4 


0.063 


1.513 


1.83 


Si2x2 


-0.901 


0.549 


1.80 


Si4x4 


-0.782 


0.668 


3.08 


OH2x2 


-1.622 


-0.173 


1.90 


OH4x4 


-0.827 


0.623 


1.90 



rec-0 2X 2 
rec-0 4X 4 



-0.432 
-0.823 



1.018 
0.625 



2.93 
3.02 



between top and bridge adsorption sites and no covalent 
binding or noticeable distortion of the graphene plane. 
Finite-size effects are clearly significant both for the bind- 
ing energies and for the chlorine equilibrium distance but 
as the energies are converging, we use the values calcu- 
lated using the 12x12 supercell for comparisons. As seen 
from the binding energies (Table [i]), chlorine physisorp- 
tion on all sites is energetically feasible in radical chlo- 
rine conditions (Eb,ci), whereas adsorption does not oc- 
cur when recombination to molecular chlorine is possible 
(Eb, c\ 2 )^ an d metastable adsorbed structures are prone 
to CI2 desorption. The lower values for the adsorption 
energy (—1.25 eV for top site adsorption calculated for 
the 12x12 supercell) in comparison to Ref. [21] are due 
to the inclusion of the additional adhesion-providing van 
der Waals correction. Using the same supercell size as in 
Ref. [21] and excluding the van der Waals correction, we 
get the chlorine adsorption energy Eb, ci = —0.76 eV/Cl, 
agreeing with the value —0.80 eV/Cl that Wehling et a/PH 
reported. 

Wu et alW^ studied one-sided graphene chlorination 
also theoretically by calculating the chlorine binding en- 
ergy on pristine graphene at different coverages. Us- 
ing atomic chlorine as the reference state, they find 
that binding is energetically stable, indicated by nega- 
tive binding energies. We test this by calculating the 
binding energy for the structure corresponding to their 
calculation, using both atomic and molecular reference 
for chlorine. For 1/4 monolayer coverage, our result 
for the binding energy with respect to atomic chlorine 
is with (without) vdW correction -0.44 (-0.29) eV/Cl, 
agreeing well with approximately —0.25 eV reported by 
Wu et al. in Fig. 4(a) of Ref. [26] Using a molecular 
reference state, however, the binding energies increase to 
1.01 (1.16) eV/Cl and the adsorbed structure is unstable, 





E B , ci =-1-867 eV/Cl 
£b,ci 2 =-0.125 eV/Cl 



E B , ci = -1-333 eV/Cl 
£b,ci 2 =0.410 eV/Cl 




E B ,c\= -1.335 eV/Cl 
E B ,c\ 2 = 0.408 eV/Cl 



Eb,c\- 
Eb,c\ 2 



-1.277 eV/Cl 
=0.466 eV/Cl 



FIG. 1. (Color online) The geometries for the different two- 
chlorine configurations, (a) CI2 molecule (b) on the opposite 
sides of the carbon layer, bridge site, bound structure [in- 
set: side view] (c) opposite sides of the carbon layer, top 
site (d) on the opposite sides of the carbon layer, bridge site, 
adsorbed structure [inset: side view]. Below the structures, 
the binding energy in radical and ambient conditions (Eb, ci 
above Eb,c\ 2 ) is given per chlorine atom. Atom colors: 
Brown / gray-chlorine, black-carbon. 



and higher in energy than our result for the adsorption 
energy of a single chlorine atom, —1.25 eV/Cl. 

More complicated adsorption structures are, of course, 
possible, as well as molecular adsorption that has pre- 
viously been studied in more detail by Rudenko et al. 
in Ref. [23] We restrict ourselves to a pair of CI atoms 
placed in the 8x8 graphene supercell. We note that if 
the distance between two chlorine atoms is initially less 
than y/7 times the carbon-carbon distance in graphene, 
they attract each other and spontaneously form chlorine 
molecules during the relaxation. Our results for the ad- 
sorption energy with respect to molecular chlorine and 
equilibrium distance of a CI2 molecule agree well with 
the results of Rudenko et al, 23 giving an equilibrium 
distance of 3.49 A and binding energy —0.250 eV/C^, 
compared to 3.58 A and —0.259 eV/Cl2 by Rudenko et 
al. in Ref. [23| The final adsorption position for a chlo- 
rine molecule is shown in Fig. TJa). Instead, for initial 
configurations where the distance between the chlorine 
atoms is longer than \fl times the carbon-carbon dis- 
tance, the mutual repulsion between the CI atoms leads 
to maximum interatomic distance between the chlorine 
atoms within the finite supercell. 

In the photochlorination experiment by Li et a/P^, 
graphene and chlorographene were deposited on the SiC>2 
substrate. In the case of graphene hydrogenation, the 
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substrate may affect the stability of adsorbed struc- 
tures, as recently reported for hydrogenated graphene, 
graphane, on SiC^P^ In order to take this possibility 
into account, we calculated the chlorine binding ener- 
gies on the a-phase Si02 surface. As Si02 is amor- 
phous, we considered four different surface termina- 
tions: oxygen-terminated (O), silicon-terminated (Si), 
hydroxyl-terminated (OH) and reconstructed oxygen ter- 
minated (rec-O) surface, for which the initial geome- 
try was taken from Ref. |32j The lattice mismatch be- 
tween the optimized substrate and graphene was only 
1.3% when a 2x2 supercell formed from graphene prim- 
itive unit cells with the optimized lattice constant cor- 
responding to freestanding graphene was used. Three 
Si02 unit cells were used in the slab, corresponding to a 
layer of 15.8 A, that was either symmetrically chlorinated 
from both sides, or terminated with an ideal graphene- 
substrate interface. The chlorine atom, graphene, and 
the three uppermost atom layers of the substrate were 
relaxed until forces acting on atoms were less than 0.01 
eV/A. The convergence criterion for the electronic de- 
grees of freedom was 10 -6 eV and a k-point mesh of 
6x6x1 was used to sample the Brillouin zone. 

For chlorine adsorption above the graphene plane, we 
find that in the radical chlorine environment, binding 
to the surface is energetically feasible on all surface ter- 
minations for the 1/8 monolayer chlorine coverage (2x2 
graphene supercell). In contrast to adsorption onto free- 
standing graphene, we find covalent bonding between the 
chlorine atom and graphene on all surfaces apart from the 
rec-0 surface, as seen in the carbon-chlorine bond lengths 
shown in Table [TTJ The absorption onto the OH surface 
is energetically most favored, as in the case of hydrogen 
adsorption.^ In the ambient environment, using molecu- 
lar chlorine as the reference state, the adsorption energy 
per chlorine atom is negative only for the OH-terminated 
surface, the adsorption energy being —0.17 eV/Cl. This, 
however, requires an ideally flat, crystalline Si02 surface 
and the effect of surface roughness is beyond the scope 
of this article. 

At lower coverages, modeled by increasing the supercell 
size to correspond to 4x4 graphene unit cells (and 2x2 
Si02 unit cells), scaling the k-point mesh correspond- 
ingly, we notice that on those surfaces (O and OH) that 
form bonds with chlorinated graphene, the absorption 
energies increase as the distance between the adsorbates 
increases (Table [H| and adsorption becomes less favor- 
able. On the rec-0 surface, on the contrary, the absorp- 
tion energy decreases. On the Si surface, graphene binds 
both to the substrate and to the chlorine atom in the 
2x2 supercell but increasing the supercell to 4x4 leads 
to change in the binding to physisorption, seen in Table |H| 
as increase in the C-Cl distance from 1.80 A to 3.08 A. As 
chlorine binds to graphene on the O and OH surfaces, the 
planarity of the carbon network is disturbed and bind- 
ing to the substrate provides additional stabilization. In 
this case, it is beneficial to have adsorbates relatively 
close to each other, keeping the carbon backbone similar 



to the chair configuration of graphane. The Si surface 
is more inert and unable to provide this stabilization at 
low coverage. On the other hand, graphene on the rec- 
O surface without chemical bonding between substrate 
atoms and the carbons resembles a freestanding graphene 
membrane with adsorbates. Like for the freestanding 
graphene, there is a repulsion between the chlorine atoms 
in neighboring periodic images, 33 and increasing the dis- 
tance between the adsorbates stabilizes the structure. At 
this lower coverage, chlorinated graphene is no longer sta- 
ble on any of the surfaces and the discrepancy between 
the experiment and theoretical calculations can not be 
simply explained by substrate effects. We exclude the 
study of chlorine adsorption on defected graphene from 
this study. 

It could, in principle, be possible that some chlorine 
intercalates between the graphene layer and the sub- 
strate through cracks and boundaries in the graphene 
layer, thus enabling chlorine attachment to graphene 
from both sides of the carbon network. Returning to 
graphene in the absence of a substrate, we probe this by 
placing two chlorine atoms on the opposite sides of the 
graphene layer, either on top of a single carbon atom, 
on neighbouring carbon atoms or symmetrically on both 
sides of a bridge site. Chlorine atoms initially placed 
on the top of neighboring carbons on both sides of the 
graphene layer relax to a bridge position. The corre- 
sponding two-chlorine geometries are shown in Figs.JTJb), 
[TJc), and[TJd) that also give the binding energy per chlo- 
rine atom in radical and ambient conditions (Eb,c\ above 
Eb,c\ 2 i respectively). The lowest-energy configuration 
corresponds to double-sided top position [ljc)] , with a 
shorter equilibrium distance than for a single adsorbed 
chlorine atom 2.5 A compared to ~ 3.1 A). The 
double-sided bridge adsorption position with a slight dis- 
tortion of the graphene layer [Fig. [ljb)] is very close in 
energy, the difference in Eb,c\ being only few meV. In- 
stead, bridge adsorption that retains graphene planarity 
[Fig. []Jd)] lies approximately 0.05 eV higher in energy. 

The binding energies per chlorine atom for the CI2 
molecules in ambient conditions (Fig. [I] below the struc- 
tures), indicate that apart from molecular adsorption, in 
ambient conditions chlorine binding in pairs on pristine 
graphene is not feasible and should be prone to desorp- 
tion processes. Being unable to explain the experimen- 
tal observations by chlorine adsorption and binding onto 
the graphene basal plane with bond formation of chlo- 
rine with sp 3 -hybridized carbon seen in the experiment, 
we thus turn our attention to the binding of chlorine to 
graphene edges. 



B. CI in zigzag graphene edges 

As a model for zigzag-terminated (ZZ) graphene edges, 
we consider an eight zigzag rows wide ZGNR (8-ZGNR, 
width ^15.6 A). For fully hydrogen-terminated ribbons, 
this corresponds to well converged edge energy with re- 
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FIG. 2. (Color online) The edge formation energy Ep,ci 2 m ambient conditions (upper panel) and binding energy per chlorine 
atom Eb,c\ (lower panel) in radical chlorine environment for the different zigzag terminations. In the upper panel, the dotted 
horizontal line shows stability with respect to pristine graphene (^b,ci 2 — 0), and in the lower panel, the red dotted line 
shows the adsorption energy of a single chlorine atom on pristine graphene in the absence of substrate (12x12 supercell, 
Ef,c\ = —1.25 eV/Cl). The geometrical structures are connected to the graph from below (above) if they are (un)stable in 
ambient conditions. Vertical dashed lines separate groups of different degree of chlorination, expressed as chlorine atoms per 
edge carbon. Brown/gray: chlorine, black: carbon, white: hydrogen. 



spect to the ribbon widthJ^l We checked this also for 
some of the chlorinated structures, finding that in the 
width range 15-28 A, the formation energies remain es- 
sentially constant. Approximately 20 A of vacuum was 
placed between adjacent ribbons to avoid interaction due 
to periodic boundary conditions. All atoms in the rib- 
bons were relaxed and the Brillouin zone was sampled 
only in the direction along the ribbon axis. Comparisons 
with systems in which the carbon atoms in the middle of 
the ribbons were fixed and only carbon atoms belonging 
to the edge hexagons were relaxed showed only energy 
differences of order 10 meV in total energies. 

We generate a large number of initial chlorinated con- 
figurations by starting both from the monohydrogenated 
(H-H-H) and the H-H-H2 edge (in which every third edge 
carbon is dihydrogenated) and by substituting some of 
the hydrogen atoms by chlorine. As all ribbon atoms 
are relaxed, some coupling between the two edges of the 
ribbon may remain. In order to take this into account, 
we also calculate all different isomers (including stereo- 
chemistry) for each edge configuration. The differences 
in total energy between the isomers are always below 0.1 
eV and in many cases they are of order meV, and thus 
the isomers are reported as a single termination, choosing 



the lowest-energy isomer as the representative structure 
for the termination. The stability of the edges is deter- 
mined by calculating the edge formation energy per unit 
length 



^Cl-xGNR 



N C E C 



graphene 



^E H2 



• NciEqi 



L 



(2) 

where £?ci-xGNR is the energy of the chlorinated ribbon, 
N a the number of a atoms, E^ raphene the energy of a sin- 
gle carbon in pristine graphene, E a the energy of species 
a and L the length of the ribbon edge. If molecular chlo- 
rine is used as the reference state, Nc\Eq\ is replaced by 
Nq\Eq\ 2 /2. A stable structure with respect to pristine 
graphene should have a negative edge formation energy. 
In order to determine whether chlorine prefers to bind 
to the edge instead of adsorbing onto the basal plane of 
graphene, we also calculate the chlorine binding energy, 



E B 



Ec 



1-xGNR 



N C E. 



ci 



graphene 



^E H2 - NciEci 



(3) 

and compare it to the lowest binding energy of a single 
chlorine atom on an infinite graphene sheet, for which 
we use the energy corresponding to the top adsorption, 
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TABLE III. The energetics for different zigzag edge termina- 
tions: the formation energy per unit length of the edge (Ef) 
and the binding energy per chlorine atom (Eb), both with 
radical chlorine (,ci) and with molecular chlorine (,ci 2 ) as the 
reference state. Also the C-Cl bond lengths are reported. The 
ordering is based on edge stability with respect to molecular 
chlorine (Ef,c\ 2 )- 



system 


Ef, ci 


Eb, ci 


Ef,c\ 2 


Eb,ci 2 


dc-c\ 




(eV/A) (eV/Cl) (eV/A) (eV/Cl) 


(A) 


H-C1-HC1 


-0.487 


-1.802 


-0.016 


-0.059 


1.88, 1.74 


H-CI-H2 


-0.245 


-1.810 


-0.009 


-0.066 


1.74 


H-H-HC1 


-0.237 


-1.755 


-0.002 


-0.012 


1.88 


H-H-CI2 


-0.458 


-1.695 


0.013 


0.049 


1.83 


H-H-H 2 


0.013 




0.013 






C1-C1-HC1 


-0.665 


-1.640 


0.042 


0.103 


1.74,1.73,1.88 


H-CI-CI2 


-0.664 


-1.638 


0.043 


0.105 


1.73, 1.83 


CI-CI-H2 


-0.425 


-1.571 


0.047 


0.173 


1.74 


H-Cl 


-0.271 


-1.335 


0.083 


0.408 


1.73 


H-H-Cl 


-0.144 


-1.066 


0.092 


0.678 


1.73 


H-C1-C1 


-0.324 


-1.200 


0.147 


0.544 


1.73 


CI-CI-CI2 


-0.728 


-1.346 


0.215 


0.397 


1.72,1.86 


CI 


-0.162 


-0.399 


0.545 


1.344 


1.75 



E B ,a = -1.25 eV. 

Table [TTT] shows the edge formation and chlorine bind- 
ing energies (^f,ci/ci 2 an ^ ^#,ci/ci 2 > respectively) of the 
chlorinated edges, along with the carbon-chlorine bond 
length (dc-ci)- Fig. [2]shows Eb,c\ m radical conditions 
(lower panel) and Ep,c\ 2 m ambient conditions (upper 
panel) , along with the geometries of structures of the ter- 
minations, depicting some of the data in Table III pic- 
torially. The structures below the horizontal lines are 
stable in their own environments (see the caption for 
details). Similarly with the lowest-energy hydrogenated 
zigzag edgesp^the lowest-energy structures are commen- 
surate with the H-H-H2 hydrogenation pattern. In all 
chlorinated edges that modify the H-H-H2 pattern, the 
binding energy per chlorine atom in radical environment 
is lower than the adsorption energy of a chlorine atom 
onto pristine graphene and thus edge binding is favored. 
Some of the chlorinated H-H-H2 modifications are stable 
also in ambient conditions (Fig. [5] upper panel), although 
the edge formation energies are very small, of the order 
-lO- 3 eV/A. The de gree of chlorination of the edge car- 
bon atoms is 2/3 in the lowest-energy configuration (H- 
C1-HC1). In general, there is only little correspondence 
between the degree of chlorination and edge stability, as 
seen in Fig. [2] 

Looking at the atomic arrangement of the lowest- 
energy structure, H-C1-HC1 (Fig. [2J bottom right), we 
notice that one of the chlorine atoms binds to a planar 



interaction of monohydrogenated ZGNRs and molecular 
chlorine and found similar stable structures with some of 
the edge carbon atoms bound to both chlorine and hy- 
drogen at low chlorine coverage of the edge atoms. The 
C-Cl bond lengths give some indication on the nature 
of binding. Although all C-Cl bond lengths agree fairly 
well with the average literature value for covalent C-Cl 
bonds of 1.76 A J™ two different regimes are to be seen in 
Table |III| The chlorine atoms are more strongly bound 



sp carbon. The carbon atom binding to both hydrogen 
and chlorine slightly relaxes upwards towards the chlo- 
rine and the chlorine is close to a top adsorption position 
for the edge carbon. Previously, Lee et a/P^ studied the 



to sp^ hybridized carbon atoms, indicated by the shorter 
bond lengths 1.72-1.75 A. The bonds between chlorine 
atoms and carbon atoms with increased sp 3 character are 
slightly longer, 1.83-1.88 A. These values are well in line 
with literature values for carbon-chlorine bond lengths in 
different chemical environments. 30 

Although we use the GNRs only to model graphene 
edges, we briefly comment on the band gaps of the chlo- 
rinated ribbons. ZGNRs with the H-H-H2 hydrogena- 
tion pattern are predicted to have a band gap that is 
also dependent on the relative arra ngement of the H2 
groups across the ribbon edges J^EU w hereas the mono- 
hydrogenated ribbons are predicted to be metallic in 
the absence of spin polarization. If spin polarization is 
taken into account, an antiferromagnetic spin structure 
between the two graphene sublattices emerges, leading 
to spin-polarized edge states in the metallic ribbonsP^ 
In the terminations commensurate with the H-H-H2 pat- 
tern, the edge carbon atoms that bind to two substituents 
are not a part of the graphene tt electron network and the 
spin-polarized edge state is destroyed. The substitution 
of edge hydrogens by chlorine atoms does not change this 
picture. The band structures in the absence of spin po- 
larization are metallic if the edge carbons are monotermi- 
nated, and in the case of a H-H-H2 modification a band 
gap is found. The presence of chlorine does not seem to 
disturb the edge state and modifies only slightly the band 
gap values. For instance, the band gap of a 8-ZGNR with 
H-H-H2 termination is 0.92 eV whereas for the H-CI-H2 
termination the band gap is 0.91 eV. 

As the realisticity of spin-polarized edge states in 
ZGNRs has been debated 37 , we note only briefly that 
the chlorinated, monosubstituted ZGNRs show spin po- 
larized edge states that are slightly lower in energy than 
the nonmagnetic ones. The spin moments on chlorine- 
binding carbon atoms are in most cases slightly enhanced 
in comparison to the hydrogen-bound ones. As an ex- 
treme example, in the H-C1-C1 terminated ribbon the 
spin magnitude at the edge carbons are 0.16 \±b and 
0.18 jiB for the carbon atoms bound to hydrogen and 
chlorine, respectively. An exception is the fully chlori- 
nated edge (named CI), where the maximal spin moment 
is only 0.06 \±b and the spin polarization is thus strongly 
suppressed in comparison to the monohydrogenated edge 
with maximal spin moments of 0.20 \±b- This may be due 
to significant chlorine contribution to the density of states 
and an extra dispersive energy band near the Fermi level 
in this structure (not shown). 
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C. CI in armchair graphene edges 

It is surprising that even though the modification of 
ZGNRs with different edge-terminating functional groups 
and the resulting effect on their electronic and mag- 
netic properties has been a field of intense study (see, 
eg., Refs. [I5j [38j and l39j) . the edge functionalization of 
armchair graphene nanoribbons (AGNRs) has caught the 
attention of the scientific community only recently!^^ 
Before that, AGNR ribbons have been assumed to be 
fully planar. Rippling of the AGNR edge in the pres- 
ence of non- hydrogen terminating groups has, h owever, 
been predicted in recent theoretical works^IlLL^ anc [ we 
take this possibility into account by using a doubled su- 
percell along the ribbon axis in our calculations. In prin- 
ciple, ripples with a longer wavelength could also be pos- 
sible, but our choice will most likely capture the essential 
physics of the system regarding the interaction between 
the edges and chlorine. Indeed, for all singly- terminated 
ribbons apart from the plain hydrogen termination, we 
found structures with edge rippling to be lower in energy 
than fully planar structures with the same termination. 

Previously, the interaction of chlorine with armchair 
edges has been only briefly studied. Lee et alW* studied 
the interaction of chlorine molecules with GNR edges but 
found only weak adsorption of chlorine molecules near 
the ribbon edge with no covalent binding. Wagner et 
a/P^ report rippling of the monochlorinated edge (the CI 
edge in our nomenclature) but they do not comment on 
its energetic stability. 

In order to study chlorinated structures of an armchair- 
terminated graphene edge, we use a 13-AGNR (width « 
14.8 A) as the model system. The electronic and ener- 
getic properties of AGNRs are known to oscillate with re- 
spect to the ribbon width ! 34 ! 40 ' Like in the case of ZGNRs, 
we consider a large number of chlorinated edge structures 
created both from the monohydrogenated and dihydro- 
genated AGNR edges, 35 and consider also different struc- 
tural isomers with respect to the two edges, as well as 
stereoisomers. Similar to the zigzag edges, the difference 
in total energy between isomers is less than 0.1 eV and 
we consider them to represent a single structure. 

The energetics of the lowest-energy chlorinated AGNR 



edges are given in Table IV , again in the order of energetic 
stability with respect to molecular chlorine (£j7ci 2 )> an d 
in Fig. [3] that also shows the geometries in the fashion of 
Fig. [2] Similar to zigzag edges, all chlorinated edges are 
stable in a radical chlorine environment in terms of edge 
formation energy Ef,cu an d in this environment, the 
chlorine binding energies are comparable for both edge 
types. In ambient conditions, however, only some of the 
AC geometries modifying the dihydrogenated edge are 
stable, and, in general, the AC edges are more stable than 
ZZ edges, having lower binding and formation energies. 
In addition, there are a larger number of stable configu- 
rations and as a matter of fact, even the doubly hydro- 
genated AC edge is more stable than the lowest-energy 
chlorinated ZZ edge. In ambient conditions, the edge 



TABLE IV. The energetics for different armchair edge ter- 
minations: formation energy per unit length of the edge (Ef) 
and the binding energy per chlorine atom (Eb), both with 
radical chlorine (ci) and with molecular chlorine (ci 2 ) as the 
reference state. Also the C-Cl bond length is reported. The 
ordering is based on edge stability with respect to molecular 
chlorine (Ef,c\ 2 )- 



system 


Ef,c\ 


Eb,c\ 


Ef,c\ 2 


Eb,c\ 2 


dc-c\ 




(eV/A) (eV/Cl) (eVA) (eV/Cl) 


(A) 




-0.990 


-2.113 


-0.173 


-0.370 


1 8^ 


HCI-H2 


-0.574 


-2.450 


-0.165 


-0.706 


1.85 


HC1 (trans) 


-0.978 


-2.088 


-0.162 


-0.345 


1.83 


H 2 


-0.091 




-0.091 






HCI-CI2 


-1.287 


-1.832 


-0.062 


-0.089 


1.80, 1.82 


H 2 -H 2 -C1 2 -C1 2 


-0.834 


-1.781 


-0.017 


-0.037 


1.78, 1.82,1.85 


H 2 -C1 2 


-0.806 


-1.721 


0.011 


0.023 


1.86 


H-Cl 


-0.342 


-1.463 


0.066 


0.281 


1.74 


Cl-H-H-Cl 


-0.327 


-1.396 


0.081 


0.348 


1.74 


H-H-C1-C1 


-0.320 


-1.368 


0.088 


0.375 


1.73 


CI 


-0.694 


-1.483 


0.122 


0.260 


1.73 



formation energy of the favored AC edge is an order of 
magnitude lower than that of the favored ZZ termination, 
-0.173 eV/A [HC1 (cis)] in comparison to -0.016 eV/A 
[H-C1-HC1]. Actually, between these two formation ener- 
gies there are four other AC edges that are energetically 
more stable than the lowest-energy zigzag edge. We do 
not find any spin polarization in chlorinated armchair 
structures. 

The degree of chlorination is higher for the most stable 
AGNR edge (1 Cl/edge carbon), compared to the most 
stable ZGNR configuration (2/3 Cl/edge carbon). Again, 
the low-energy structures modify the H2 edge and have 
sp 3 hybridized carbon atoms with the C-Cl bond lengths 
1.83-1.85 A. The bonds between sp 2 carbons and chlorine 
are slightly stronger and thus shorter, the bond length 
being approximately 1.73 A. Notably, the binding energy 
per chlorine atom in ambient conditions for the lowest- 
energy structures (Table |lV| Eb,c\ 2 ) is n °t only lower 
than chlorine atom adsorption (0.499 eV/Cl) but also 
lower than the adsorption energy of molecular chlorine 
(-0.125 eV/Cl). 



D. Thermodynamic stability 

The DFT-calculated energies do not take into account 
the environment, for instance the pressure of the gaseous 
species and the finite temperature. It is intuitively clear 
that if there is only a very limited amount of a reagent 
available, this may have an effect on the favored prod- 
uct. We may take t hese effe cts into account using ah 
initio thermodynamics^ 2 * 35 1 411 that allows us to compute 
the relative stabilities of different edge terminations as a 
function of the chemical potential of the gaseous species, 




T<?m TM TOT 



FIG. 3. (Color online) The edge formation energy Ep,ci 2 i n ambient conditions (upper panel) and binding energy per chlorine 
atom Eb, ci (lower panel) in radical chlorine conditions for the different armchair terminations. In the upper panel, the dotted 
horizontal line shows stability with respect to pristine graphene (Eb,c\ 2 — 0), and in the lower panel, the red dotted line 
shows the adsorption energy of a single chlorine atom on pristine graphene in the absence of substrate (12x12 supercell, 
Ef,ci = —1.25 eV/Cl). The geometrical structures are connected to the graph from below (above) if they are (un) stable in 
ambient conditions. Vertical dashed lines separate groups of different degree of chlorination, expressed as chlorine atoms per 
edge carbon. Atom colors as in Fig. [2] 



the chlorine and the hydrogen. Although gaseous hy- 
drogen is most likely scarcely available in the reaction 
conditions the resulting phase diagram can be used as 
a qualitative rough guide to the effects of other forms of 
hydrogen, for instance hydrogen bound to the substrate. 
The Gibbs energy of the edge formation is given by 



AG = E FC \ 2 



(4) 



where the reference state for chlorine is chosen to be 
molecular chlorine, and pi gives the edge density of 
species i and \±i the chemical potential of species i. By 
using the formation energies calculated with respect to 
atomic chlorine and substituting pciMci 2 /2 by pciMci, the 
Gibbs energy with respect to atomic chlorine can be de- 
termined. The chemical potential, given by 

Mi (T) = H°(T)-H°(0)-TS°(T) + k B Tln(P i /P°), (5) 

where H° (S°) is the enthalpy (entropy) in standard pres- 
sure and Pi the partial pressure of species z, can b e can 
be calculated using thermodynamical reference data.^E^ 
Fig.[4]shows the chemical potential of gaseous atomic and 
molecular chlorine calculated using Eq. ([5| at different 
temperatures and partial pressures. 
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FIG. 4. (Color online) The chemical potential of molecu- 
lar (violet/dark gray) and atomic (green/light gray) chlorine 
(a) as a function of temperature at three values of P (b) as a 
function of partial pressure at room temperature (T=298 K). 
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FIG. 5. (Color online) The phase stability diagrams showing the thermodynamically favored edge termination (with lowest 
AG) in different hydrogen and chlorine environments, given by the chemical potential of the gaseous species CI2, CI and H2. 
(a) Zigzag edges (b) Armchair edges (c) All edges (d) All edges including pristine graphene. 



As noted in the previous sections, in general the 
armchair-type terminations are more stable than their 
zigzag-type counterparts. Figs. [5|a) and|5|b) show sep- 
arately the stability diagrams of the zigzag and armchair 
edges, that is the edge termination with the lowest AG, 
as a function of the hydrogen and atomic chlorine chem- 
ical potentials. For comparison, also the chemical poten- 
tial scale corresponding to molecular chlorine is given. 
As the chemical potential, and thus the partial pressure, 
of chlorine increases, the more densely chlorinated edge 
terminations become favored. Fig. [5|c) combines subfig- 
ures (a) and (b). Zigzag edges are thermodynamically 
the most stable structures only at low chlorine and hy- 
drogen chemical potentials. 

If AG of a termination is negative, an infinite graphene 
layer can be spontaneously broken into edgesP^MI Thus, 
the possibility of pristine graphene being the favoured 
structure should be taken into account in the stability 
diagrams. Fig. j5^d) shows a combined armchair-zigzag 
stability diagram including a region for pristine graphene. 
At low hydrogen and chlorine chemical potentials, that is 
at low pressures or concentrations, edge formation is not 
possible and pristine graphene is preferred. All chlorine- 
or hydrogen-containing terminations that are stable with 
respect to pristine graphene are of armchair type, and 
they are only favored over graphene at higher hydrogen 
and chlorine chemical potentials. 

Our approach does not consider chemical kinetics, or 
reaction barriers between different structures. In a real 
environment, metastable states are possible and this 
might allow some of the seemingly unstable chlorinated 
structures also to exist in ambient conditions. It is dif- 
ficult to asses the experimental conditions used by Li et 
a/PO as they do not estimate their chlorine partial pres- 
sure. They, however, mention that the reaction is per- 
formed at room temperature, at 1 atm pressure in a mix- 



ture of radiated CI2 and N2 so this gives an upper limit 
to the chlorine partial pressure. As seen in Fig. |4jb), at 
room temperature this corresponds to a maximum chem- 
ical potential of approximately —0.45 eV that would give 
the armchair HCI-CI2 as the most stable edge configura- 
tion and allow its formation from pristine graphene. 



IV. DISCUSSION AND CONCLUSIONS 

We have studied the interaction of chlorine with 
graphene using density-functional theory. By comparing 
the energetics of different chlorine-containing edge ter- 
minations, both for armchair and zigzag graphene edges, 
to chlorine adsorption on pristine, freestanding graphene 
as well as graphene on the Si02 surface, we showed that 
chlorine binds to the edges rather than adsorbs onto the 
basal plane. In general, armchair edges are more read- 
ily chlorinated and the resulting structures are lower in 
energy compared with both atomic and molecular chlo- 
rine adsorption on pristine graphene. In addition, using 
ab initio thermodynamics, we studied the effect of en- 
vironment chlorine and hydrogen content on the stabil- 
ity of different edges and found that pristine graphene 
may be spontaneously broken into edges in a sufficiently 
chlorine-rich environment. The presence of the silicon 
dioxide substrate lowers the adsorption energies onto the 
basal plane of graphene but even then the binding to the 
edges is energetically favored. 

In the experiment, Li et a/P^ suggest that chlorine 
binds covalently to the basal plane of graphene, using 
the increase in the disorder-induced Raman D-peak as 
evidence, as well as XPS spectra showing covalent bonds 
between sp 3 hybridized carbon atoms and chlorine atoms. 
We suggest that the nanodomains observed in the exper- 
iment by Li et alW* may actually be fractured graphene 
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with chlorine-terminated edges. Our interpretation by 
no means contradicts these experimental observations. 
An enhanced reactivity upon chemical modification at 
graphene edges has previously been reported.^ In addi- 
tion, the Raman D-band is activated also through scat- 
tering from graphene edges Carbon atoms that are 
sp 3 -hybridized and bind to chlorine are present in the 
lowest-energy edge terminations. As seen in our results in 
the presence of a substrate, chlorine attachment onto the 
basal plane may be locally feasible if additional stabiliza- 
tion through the substrate is available. Local substrate- 
stabilized covalent binding could also explain the cover- 
age, 8 atom-%, estimated by Li et al based on a XPS 
measurement that seems quite high to be achieved only 
through chlorinated edges. 

Li et alW^ also found that the average height of the 
chlorinated graphene increases from 0.9 nm to 1.1-1.7 nm 
upon the reaction with chlorine. As the carbon-chlorine 
bond length is only approximately 0.18-0.35 nm, the ex- 
tremes corresponding to covalent bonds and physisorp- 
tion, binding to chlorine alone can not explain the in- 
creased height and roughness. It is possible that the for- 
mation of nanodomains is also related to surface mor- 
phology, so that graphene on higher regions or in re- 
gions with a certain surface termination react more read- 
ily, thus increasing the height variation. Additionally, 
in fractured graphene, chlorine may be able to interca- 
late between the substrate and graphene, increasing the 



substrate-carbon distance. Yet another possibility is that 
during the edge formation, the interaction with chlorine 
leads to upward bending of the edge regions, or partial 
overlap between adjacent edges. Cleavage only in the 
high-symmetry armchair and zigzag directions is unlikely 
to occur in the reality, and also bulky chlorinated hy- 
drocarbon groups [of type (CH x Cl2- x )nCH y Cl3_ y ] could 
form at the edges. Even if the Raman maps of the 
chlorine-reacted graphene do not show spatial variation, 
the resolution of the measurement with a beam of 1 /im 
in diameter is much larger than the lateral dimension of 
the observed domains, 30-50 nmP 

Our results show that alternative explanations than 
the simple basal plane adsorption, such as a combination 
of chlorine-containing graphene edges and substrate ef- 
fects studied in this work are needed to explain graphene 
chlorination in experiments. Additionally, defects are 
likely to be important as reactive sites in the basal plane 
but this is a topic of future work. 
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